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In  this  work  we  study  the  annihilation  of  dark  matter,  considering  it  as  a  neutral  particle  with  magnetic 
and/or  electric  moments  not  null.  The  calculation  of  the  effective  section  of  the  process  XX  >  77  is 
made  starting  from  a  general  form  of  coupling  XX7  in  the  framework  of  an  extension  of  the  Standard 
Model.  We  found,  when  taking  into  account  an  annihilation  of  DDM-antiDDM  to  monoenergetic 
photons,  that  for  small  masses,  m  <0  GeV,  an  electric  dipole  moment  ~  10-6  e  cm  is  required  to  satisfy 
the  current  residual  density,  while  for  the  range  of  greater  sensitivity  of  HAWC,  10  TeV  <  E  <  20  TeV, 
the  electrical  dipole  moment  must  be  of  the  order  of  10  8  e  cm. 
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1.  Introduction 

The  enigma  of  dark  matter  is  perhaps  the  most  interesting 
problem  of  modern  astrophysics,  so  much  so  that  it  has  led 
to  the  incursion  of  elementary  particle  physics.  The  joint 
work  of  these  two  disciplines  has  as  one  of  its  main  objectives 
to  determine  the  nature  and  properties  of  dark  matter, 
either  through  direct  or  indirect  detection.  This  enigma  of 
the  missing  mass  has  been  a  problem  since  Zwicky  in  1 933 
measured  the  masses  of  extragalactic  systems  [1].  Nowadays, 
given  the  evidences  of  the  galactic  dynamics  (rotation  curves), 
galaxy  clusters,  structure  formation,  as  well  as  the  Big  Bang’s 
nucleosynthesis  and  the  cosmic  background  radiation, 
it  is  suggested  that  baryons  can  only  explain  matter,  the 
majority  of  the  missing  mass  must  be  non-baryonic.  The 
non-baryonic  nature  of  dark  matter  is  clear  evidence  that 
our  understanding  of  the  matter  components  of  elementary 
particle  physics,  beautifully  described  by  the  Standard 
Model  (SM)  is  incomplete.  For  this  reason,  theoretical 
physicists  have  considered  new  physics  beyond  the  Standard 
Model  in  order  to  accommodate  (at  least)  a  non-baryonic 
candidate  as  dark  matter  (DM),  since  the  only  dark  matter 
candidate  in  the  SM  is  the  neutrino,  which  is  inadequate 
to  explain  most  of  the  DM  [2].  The  most  promising 
candidates  that  emerge  beyond  the  Standard  Model,  are  the 
massive  particles  of  weak  interaction,  commonly  known  as 


WIMPs,  examples  of  these  are  the  neutralino  [3]  and  the 
axion  [4],  but  unfortunately  they  have  not  been  detected. 
In  the  absence  of  the  discovery  of  such  particles  it  is  worth 
exploring  other  possibilities.  An  alternative  line  of  research 
is  to  take  an  approach  independent  of  the  model  and  try 
to  phenomenologically  explore  the  possible  properties 
of  a  dark  matter  particle.  On  this  line,  the  restrictions 
for  strongly  interacting  dark  matter  were  considered  in 
Ref.  [5].  In  addition,  the  autointeraction  of  dark  matter 
has  been  considered  in  the  Refs.  [6,7].  Some  people  have 
studied  whether  dark  matter  could  be  charged  [8]  or  have  a 
milichargue  [9, 1 0] .  Likewise,  it  has  been  studied  that  within 
these  phenomenological  possibilities,  dark  matter  has  an 
electric  or  magnetic  dipole  moment  [11,  12,  13-13].  In  this 
work  we  consider,  precisely,  the  possibility  that  dark  matter 
possesses  an  electrical  and/or  magnetic  dipole  moment  and 
can  emit  gamma  radiation  as  a  result  of  its  annihilation. 
In  the  next  section,  we  introduce  the  effective  Lagrangian 
for  the  interaction  of  dipolar  dark  matter  with  photons. 
The  calculation  of  the  annihilation  cross-section  by  the 
thermally  averaged  relative  velocity  is  calculated  in  section 
3.  In  section  4  we  speak  of  monoenergetic  gamma  radiation, 
while  in  section  5  conditions  are  established  on  the  magnetic 
and  electric  dipole  moment  to  satisfy  the  current  residual 
abundance.  We  give  our  conclusions  in  section  6. 
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2.  Dipolar  Dark  Matter 

Dark  matter  seems  to  be  made  up  of  non-relativistic  particles 
that  interact  only  gravitationally  and  perhaps  by  weak 
interaction.  In  general,  the  coupling  to  photons  is  assumed  to  be 
non-existent  or  very  weak,  so  that  electromagnetic  interactions 
have  not  been  seriously  considered  [12].  However,  although 
the  DM  particles  are  considered  chargeless,  that  is,  they  have 
no  electrical  charge,  they  could  be  coupled  to  photons  through 
loops  in  the  electric  and  magnetic  dipole  moment  [13].  For 
this  reason,  we  assume  that  dipolar  dark  matter  is  a  class  of 
WIMPs,  which  are  endowed  with  dipole  moments,  so  they  can 
interact  weakly  via  electromagnetic  interaction.  The  particles 
proposed  as  DM  are  Dirac  fermions,  that  is,  spin  particles 
1/2,  because  for  a  permanent  electric  and/or  magnetic  dipole 
moment  the  particle  must  have  a  spin  other  than  zero  [15]. 
The  phenomenology  of  the  dark  matter  particles  charged  is 
determined  to  a  large  extent  by  the  ability  of  these  particles 
to  form  bound  states,  such  as  atoms,  with  electrons,  protons 
or  each  other.  Interestingly,  a  neutral  particle  with  an  electric 
dipole  moment  (D)  can  make  a  state  linked  with  an  electron, 
as  observed  by  Fermi  and  Teller  [16],  but  only  if  the  dipole 
moment  is  greater  than  0.639e  a0  =  3.4x  10  9e  cm  (assuming 
that  m  »  m ),  where  an  is  the  radius  of  Bohr.  For  small 

g  er  0 

dipole  values,  the  electron  sees  both  poles  of  the  dipole  and  is 
in  an  unstable  orbit.  This  critical  electric  dipole  moment  scales 
inversely  with  the  reduced  electron-dipole  mass,  so  a  state 
bound  to  a  proton  can  occur  if  the  mass  of  the  dipole  is  »m^ 
and  D  »  1.8  X  10“12  e  cm,  however,  such  D  values  can  not 
occur  for  point  dipole  dark  matter.  Also,  the  weakness  of  the 
dipole-dipole  interaction  prevents  the  formation  of  stable  dark 
matter  atoms  [15].  On  the  other  hand,  if  the  particles  of  dipolar 
dark  matter  (DDM)  are  the  dark  matter  of  our  Universe,  had 
to  undergo  a  thermal  history  similar  to  that  of  any  WIMP,  that 
is,  existed  in  the  early  Universe  when  the  temperature  T  >  m  , 
and  their  interactions  froze  when  T  fell  below  resulting  in 
some  residual  cosmological  abundance  [17].  For  a  weak-scale 
thermal  particle,  the  residual  abundance  in  the  case  of  s  wave 
annihilation  is  approximately  established  by  [18] 


Qh2 


3x10  26  cm3^s 

Gnn*') 


(1) 


examples  of  these  experiments  are  CRESST,  XENON, 
CDMS,  DAMA  and  COGENT.  Instead  indirect  methods 
allow  us  to  detect  the  WIMP  through  the  observation  of 
the  products  emitted  in  their  annihilation  in  the  galactic 
halo  or  in  the  interior  of  the  Sun  and  the  Earth  where  they 
could  have  been  gravitationally  trapped.  In  this  annihilation 
certain  types  of  radiation  would  be  emitted,  such  as:  high 
energy  photons  (gamma  rays),  neutrinos,  electron-positron 
pairs,  proton-antiproton  pairs,  among  others.  An  example 
of  this  type  of  experiments  is  HAWC  (High  Altitude  Water 
Cherenkov).  This  gamma  ray  observatory  is  built  on  a  site 
located  4, 1 00  meters  above  sea  level  on  the  northern  slope 
of  the  Sierra  Negra  volcano,  which  is  located  in  the  State  of 
Puebla,  Mexico.  It  is  a  set  of  300  Cherenkov  water  detectors 
sensitive  to  gamma  rays  and  charged  particles  with  energies 
from  100  GeV  to  100  TeV  [19]. 

2. 1  The  Effective  Lagrangian  for  Coupling 

The  effective  Lagrangian  for  the  coupling  of  a  Dirac  fermion 
with  a  magnetic  dipole  moment  and  an  electric  dipole 
moment  with  the  electromagnetic  field  is 

+  (2) 

For  low  energies  compared  to  the  mass  of  dark  matter,  the 
photon  is  blind  to  the  differences  between  M  and  D.  For  the 
interaction  of  Equation  (2),  pairs  of  DDM-antiDDM  particles 
in  the  galactic  halo  or  contained  in  any  region  of  the  Universe 
with  high  densities  (centers  of  galaxies,  clusters  of  galaxies),  can 
annihilate  directly  to  gX,  where  X  can  be  some  other  neutral 
state,  such  as:  7,  2T  or  a  Higgs  boson.  In  this  work  we  assume 
that  the  annihilation  of  DDM  particles  is  towards  two  photons 
through  the  diagrams  shown  in  Figure  1.  Said  annihilations 
take  place  through  s  waves,  so  (J  v  is  almost  independent  of 
the  speed  and  therefore  independent  of  the  temperature  [17]. 


Said  mass  density  of  residual  DDM  particles  is  fixed  by  the 
annihilation  cross  section,  cr  ,  of  all  the  lighter  particles 
multiplied  by  the  relative  velocity  v.  Note  that  smaller 
effective  annihilation  sections  correspond  to  much  higher 
residual  densities.  For  the  detection  of  DDM  particles 
(WIMPs),  supposedly  integrating  the  galactic  halo,  both 
direct  and  indirect  methods  are  used.  To  detect  the  WIMP 
by  the  direct  method  is  done  by  the  nuclear  recoil  that 
they  produce  in  their  elastic  dispersion  with  the  nuclei  of 
the  detector  used  as  target  in  the  laboratory  [11,  12,  13], 


Figure  1.  Feynman  diagrams  for  the  annihilation  of  a  pair  of 
DDM-antiDDM  to  two  photons. 

3.  Calculation  of  the  Effective  Section  of  the 
Annihilation  Process:  XX  D 

We  consider  the  process  of  annihilation,  at  tree  level,  of 
dipolar  dark  matter  (DDM)  to  two  photons.  The  calculation 
of  the  cross  section  is  made  from  the  frame  of  reference  mass 
center  (MC).  It  is  proposed  as  a  propagator  to  the  dipolar 


ISSN  No.:  2321-8649(Print)  ISSN  No.:  2321-9289(Online);  Registration  No. :  CHAENG/20 18/5 1628 


J.  Nucl.  Phys.  Mat.  Sci.  Rad.  A.  Vol.  6,  No.  1,  Aug.  2018 


pp.35 


dark  matter  itself.  For  this  process  we  have  two  contributions 
of  low  order  (see  Figure  2),  so  the  amplitude  of  annihilation 
of  dipolar  dark  matter  will  be  M  =  Mx  +  Mr  The  amplitude 
of  the  first  diagram  is  given  as 


Mx  =  -*2„  V;(4&(*i)[u  (p2)avp(M  +  ^D)\ 

xia^iM  +  ^D)  ufo)], 


^1+Wx 


while,  the  amplitude  of  the  second  diagram  is 


So,  the  total  amplitude  M  considering  the  two  diagrams  for 
this  process  is 
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Figure  2.  diagram  1  at  the  top  and  diagram  2  at  the  bottom. 


where  £lm(v)  and  £]y(m)  are  the  components  of  the  four- 
moment  kx  and  respectively;  qx  and  q2  are  the  transferred 
moments  and  m  is  the  mass  of  the  particle  of  dipolar  dark 
matter.  Before  obtaining  the  square  of  the  amplitude  \M  \, 
some  simplifications  can  be  made,  since  p2x  =  ^  =  0, 

and  IP2  =  0  given  that  the  resulting  particles  are  photons,  the 
denominators  of  the  propagators  are  q2  —  mx  =  —  2px  •  kx 

and  ql —m2x= -2px- k2.  In  addition,  Ward’s  identity 
is  used,  ^2  u( Pi  X  u( P\  X  =  ( A  +  mx  Xc  >  as  weH  as  the 
relationships  of  completes  ^2 p  ie*n£v  ~ an<i 

u(  y>2  u(  p2)a  =  (^2  —  mx  )da  .  Considering  all  the  above 
we  have  that  the  square  amplitude  averaged  over  the  spins  is 


Wf=\Tr 


gm^h~mx 


klvk^ 


Pi  'K 


ov,>  (M  +  75  D)(^  +mx)a 


n\ 


x(M  +  75Z))  +  - 


'  1^  2  n  u\ 
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Pi-K  Pi  ’4 


(6) 


The  calculation  of  this  last  trace  results  in  199  terms,  but 
when  taking  into  account  that  the  particles  resulting  from 
this  process  are  photons,  then  the  terms  that  have  a  factor  a 
IP  and  IP  will  be  canceled  because  IP  =  IP2  =  0.  Therefore,  the 
trace  is  reduced  to  only  33  terms  expressed  as  four-moment 


products  p  ,  p2 ,  kx  and  k  ,  which  at  the  same  time  we  can 
write  in  terms  of  the  basic  kinematic  variables:  energies  and 
dispersion  angle  of  the  CM  system.  Finally,  in  terms  of  these 
variables  and  with  a  bit  of  algebra,  we  arrive  at  the  following 
expression  of  the  square  amplitude  averaged  over  the  spins: 


\m\2  = 


n2  2  2  n 

h  —  p  cos  6 


16 


-mx  +  p2m2  +  Ex  +m2£2  —  (2p 2m2  +  2p2£2)cos2  6 

i  4  Aq 

+p  COS  u 
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!  +  3p2^2  -Ex  +  3m2E2  +(— 6p2 


+2p2^2 )  cos2  0  —  p4  cos4  6 


M2D2 


(M4  +  £>4) 


(7) 


where  E^  and  m  ,  are  the  energy  and  mass  of  the  dipolar 
dark  matter  (or  antiDDM),  respectively,  p  is  the  moment 
of  the  incident  dark  matter  particles  and  0  is  the  scattering 
angle  of  the  outgoing  particles. 


Because  there  are  only  two  particles  in  the  final  state 
and  we  are  working  in  the  frame  of  reference  of  the  center  of 
mass,  we  will  use  the  Golden  Rule  of  dispersion  for  the  case 

1  +  2  ->  3  +  4: 
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dcr 

dfl 
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64t r2  (£,  +E2)2  |Pi 
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For  our  problem  and  based  on  the  conservation 
of  momentum  and  energy,  we  have  to:  I  Pi  |=  P , 


Ei  =E2  =Ex  =  ,/p2+ot2  ,  Ei=Ei=\v3\  =  Ex, 

h  -  v2 1 1=  vrd  =  2p / Ex  ,  with  Ex  =  (|  p,  |2  +«22  )‘/2  . 
Substituting  the  above  in  Equation  (10),  we  find  that  the 
differential  effective  section  of  annihilation  of  dipolar  dark 
matter  to  photons  is: 


d<JM 
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64tt2  (2E  YpE2- p2  cos2  0 
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+  32(ra4  +  3p2^2  -E4  +3m2E2  +(-6p2m2  +2p2E2)cos2  6- 


4  4 

p  COS 


o)m2d2 


(9) 


Then,  taking  into  account  that  the  process  is  symmetric  with 
respect  to  the  collision  axis,  the  phase  space  is  written  as  an 
integral  over  the  polar  angle  6  of  the  reference  frame  of  the 


CM,  so  when  performing  the  corresponding  integrals  we 
have  to: 


=  - 


67T 


[m4  +d4) 


Gy2  Arc  Tan  {^[2  — 1  j 

A2-1 
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+  M2D2 


12  y2  Arc  TaniyJ  y2  —  1  j  ^ 

- - —  T  34 
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where  y  =  mx  j  Ex  . 

Finally,  if  we  want  to  know  about  the  average  of  the 
thermal  distribution  of  the  WIMPs  we  do  it  by  means  of  the 
effective  section  of  annihilation  by  the  thermally  averaged 
relative  velocity,  (crannvreY)  .  Because  we  are  working  within 
the  non-relativistic  limit,  it  is  possible  to  use  the  method 
described  in  Ref  [20],  which  qualitatively  allows  us  to 
perform  the  integrals  to  obtain  a  ann vrd  and  then  expand 
on  vreP  we  expand  on  vn[  and  then  we  make  the  remaining 
integrals  on  analytically  simple  trigonometric  functions  to 
arrive  at  the  final  answer: 

Ka/>  =  1.7l423xl034cm3s-1 

x  [6  ( AT 4  +  GM2D2  +  D 4 ) 

+  (3Af 4  +  2M1  D1  +  3D4 )  (v2el )]  (11) 

So  that,  using  (<Jannvrel )  &  ard  +  brd  (v2rtl  )  =  ard+  Gbrd  /x 

[21]: 

((j^^  )  =  i  .71423  x  1 034  cm3s_1 

6{M4  +6M2D2  +D4) 
x^Gevx  +(3Af4  +2M2D2  +3T>4)  — 

(12) 

where  mGeVx  =  / GeV  y  x  =  mx  /T  is  the  decoupling 
energy  that  within  this  non-relativistic  limit  x  »  1  (o  T 
«  for  WIMPs  x  =  22  [22]. 


4.  Monoenergetic  Gamma  Radiation 

The  DDM-antiDDM  pairs  can  be  annihilated  to  two 
photons  through  the  diagrams  shown  in  Figure  1.  Because 
the  particles  move  with  speeds  below  one  thousandth  part 
the  speed  of  light,  v  =  230km  s-1 ,  the  photons  produced 
will  be  very  closely  monoenergetic  with  energies  equal  to  the 
resting  mass  of  the  DDM  particles  E  =  mx  ,  which  allows 
the  gamma  rays,  in  the  GeV  and  TeV  range,  to  access  the 
most  relevant  mass  range  of  dark  matter  particles.  Gamma 
rays  are  ideal  for  studying  the  annihilation  of  dark  matter, 
because  they  do  not  deviate  in  intermediate  magnetic  fields 
and,  therefore,  point  to  the  place  where  they  are  created. 
This  allows  us  to  look  for  gamma-ray  signatures  not  only 
in  our  neighborhood  of  the  galaxy,  but  also  in  distant 
objects  such  as  satellite  galaxies  or  even  clusters  of  galaxies. 
Another  advantage  of  the  use  of  gamma  rays  is  that,  in  the 
local  Universe,  gamma  rays  do  not  suffer  attenuation  and, 
therefore,  retain  the  source  spectral  information  intact  on 
Earth  [23].  These  characteristics  of  gamma  rays  make  the 
HAWC  observatory  ideal  for  studying  candidates  for  dark 
matter,  specifically,  dipolar  dark  matter.  Although  FiAWC 
is  sensitive  to  photons  from  100  GeV  to  100  TeV,  it  has  a 
maximum  sensitivity  in  the  range  of  10  to  20  TeV,  which 
makes  it  sensitive  to  diverse  searches  for  dark  matter 
annihilation,  including  extended  sources,  emission  diffuse 
of  gamma  rays,  and  the  emission  of  gamma  rays  of  subhalos 
of  non-luminous  dark  matter.  A  subset  of  these  sources 
includes  dwarf  galaxies,  galaxy  M31,  the  Virgo  cluster 
and  the  galactic  center.  Likewise,  the  response  of  FiAWC 
to  gamma  rays  from  these  sources  has  been  simulated  in 
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several  channels  of  well-motivated  dark  matter  annihilation  an  annihilation  of  DDM-anti  DDM  to  monoenergetic 

(bb,  tt,r+r~ ,W+W~)  [24].  photons  (see Table  1). 


5.  Results 

For  our  analysis,  we  consider  that  D  =  M,  so  that  the 
expression  of  the  effective  section  of  annihilation  by  the 
relative  thermally  averaged  speed  for  the  process  XX  >  77 
is: 


(aannvrel  }  =  1 .71423  x  1034  cm3s 


3-12 
"  m. 


'GeVX 


48  Da 


48Z)4 


(13) 


We  are  taking  into  account  the  upper  limit  for  dipolar 
moments,  D  =  M  <  3  x  10-16  (e  cm),  reported  by  K. 
Sirgurdson  et  al.  (Ref.  [13])  and  the  decoupling  energy  for 
the  WIMPs,  x  =  22  ,  we  can  write  the  previous  expression 
as  follows 


(Vanned  )  =  0-66649  X  10  26  CHl’s 


D 

A«+— 

22 


(14) 


where  D*6  =  (Z)  /3  x  10-16)4 .  Note  that  ( <Jannvrei )  has 
the  order  of  magnitude  corresponding  to  the  total  self- 
annihilation  cross  section  for  a  generic  WIMP,  which  is 
usually  set  as  (a  v)  «  3  x  10_26cm3s_1  [25]. 

On  the  other  hand,  since  the  residual  density  of  cold 
dark  matter  estimated  by  WMAP  is  Vih1  «  0.11  [26],  we 
can  relate  Equation  (14)  with  Equation  (1)  to  obtain: 


Table  1:  Required  masses  of  dipolar  dark  matter  and  range  of 
values  of  electrical  dipole  moment  that  satisfies  residual  abundance 

Vlh2  =0.11. 


Mass  range  m  [GeV] 

Range  of  D  [e  cm] 

6  <  m  <10 
—  x  — 

2.37  x  10-16  <  D  <  3  x  10-6 

10  <mx<  100 

7.5  x  10"17  <  D  <  2.37  x  Vy* 

100  <  mx<  1000 

2.38  x  10-17  <  D  <  7.5  x  10"17 

1000  <m  <  100000 
—  x  — 

2.38  x  10-18  <D<  2.37  x  lO'17 

10  TeV  <mx<  20  TeV 

5.3  x  10"18  <  D  <  7.5  x  1018 

0.20 

0.15 

\  010 

0.05 

0.00 

Figure  3.  The  graph  shows  the  residual  abundance,  Equation  (13), 
for  D  =  3  x  10-16  e  cm 


Qh2 


3x10  26  cm 3/s 

K  A/) 


0.11. 


(15) 


Then,  if  DDM  particles  are  cold  dark  matter,  Equation  (15) 
implies  that  the  DDM  particle  mass  density  depends  on 
the  electrical  dipole  moment  D  and  the  mass  therefore, 
for  electrical  dipole  moments  below  <  3  X  10“16  e  cm  and 
masses  >  6  GeV  it  is  always  satisfied  with  the  residual 
abundance  given  by  WMAP  (see  Figure  3). 

The  value  of  the  electrical  dipole  moment  that  causes 
equality  to  be  met,  in  the  mass  range  of  the  DMM  particles 
6  GeV  <  mx  <100  TeV  ,  becomes  up  to  126  times  more 
smaller  than  D  =  3  X  10“16  e  cm,  as  can  be  seen  in  Figure 
4.  But,  for  values  in  the  small  mass  range,  GeV,  an  electric 
dipole  moment  of  the  order  of  ^T0“18  e  cm,  is  required 
to  satisfy  the  current  residual  abundance,  in  accordance 
with  Ref.  [8];  while  in  the  range  of  greater  sensitivity  of 
HAWC,  1 0  TeV  <  m  <  20  TeV,  the  electric  dipole  moment 
must  be  of  the  order  of  ^10“18  e  cm,  taking  into  account 


Figure  4.  The  graph  shows  the  range  of  values  of  both  D  and  , 
of  in  which  the  residual  abundance  is  satisfied  Qh2  =  0.11 . 
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6.  Discussion  and  Conclusions 

In  this  work  we  study  the  annihilation  of  dark  matter, 
considering  it  as  a  neutral  particle  with  magnetic  and/or 
electric  moments  not  null.  For  this  reason,  a  neutral  WIMP 
with  an  associated  magnetic  and/or  electrical  dipole  moment 
of  the  order  of  D  =  M  <  3  x  10“16  e  cm  is  proposed.  The 
effective  section  of  annihilation  XX  ~ ^  77  was  calculated  in 
analytical  form.  Likewise,  this  effective  section  was  analyzed 
numerically  (< &annvrd },  where  it  was  found  that  for  small 
masses,  <  0  GeV,  an  electrical  dipole  moment  M0“16  e  cm, 
is  required  to  satisfy  the  current  residual  abundance,  while 
considering  that  HAWC  is  sensitive  to  the  photons  of  100 
GeV  at  100  TeV  we  need  a  dipole  electric  moment  D  <  7.5 
X  10"17  e  cm.  In  addition,  it  was  found  that  in  the  HAWC’s 
highest  sensitivity  range,  10  TeV  <  <  20  TeV,  the 

electrical  dipole  moment  must  be  of  the  order  of  10“18  e  cm. 
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